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We calculate the corrections to the Higgs mass in general theories restricted to the case of mass¬ 
less gauge bosons (the gaugeless limit). We present analytic expressions for the two-loop tadpole 
diagrams, and corresponding expressions for the zero-momentum limit of the Higgs self energies, 
equivalent to the second derivative of the two-loop effective potential. We describe the implemen¬ 
tation in SARAH, which allows an efficient, accurate and rapid evaluation for generic theories. In the 
appendix, we provide the expressions for tadpole diagrams in the case of massive gauge bosons. 


I. INTRODUCTION 

With the discovery of the Higgs boson in the range of 125 126 GeV the standard model (SM) has been completed 

PJ[2]. The uncertainty in the Higgs mass measurement has continuously decreased and is well below 0.5 GeV today [3]. 
This small uncertainty is currently much better than the theoretical prediction in any scenario beyond the Standard 
Model (BSM). Therefore, more precise calculations are necessary to better confront BSM models with the Higgs mass 
measurement. This has two motivations of particular weight: (1) In the minimal supersymmetric standard model 
(MSSM) and many extensions thereof, radiative corrections are required to be at least as significant as the tree-level 
contribution, so higher-order corrections are especially important. (2) In the Standard Model and non-supersymmetric 
extensions thereof a precise calculation is required to extract the parameters of the model, which when run to high 
energies gives information about the stability or lifetime of the potential - which may point the way to new physics if, 
as appears currently, the potential is metastable. Beyond these motivations, for a generic model of new physics with 
boundary conditions fixed from the top down (such as supersymmetric models) it is important to know what regions 
of parameter space are allowed, compatible with the measured Higgs mass. For example, a one-loop calculation may 
naively lead to excluding certain constrained scenarios, whereas with a two-loop calculation the Higgs mass may be 
large enough; this is related to the difficulty in estimating the error in the Higgs mass calculation, since at two-loop 
order there are new contributions from particles that have no direct couplings to the Higgs, and a simple variation of 
the renormalisation scale as an estimate of the error is not sufficient. 

In general, there are three approaches to tackle the problem of finding the Higgs mass precisely: (i) effective potential 
calculations, (ii) diagrammatic calculations, (iii) renormalisation group equation methods. We shall concentrate in 
the following on the first two options. Calculations from the effective potential suffer from a larger uncertainty 
compared to diagrammatic calculations because of the missing momentum contributions. However, these are only 
really pronounced at one-loop level, and it is already possible to calculate the full one-loop Higgs mass inlcuding 
momentum dependence for generic models using SARAH nn to produce SPheno B ITO] output or SOFTSUSY EH 
output via FlexibleSUSY [12]; explicit results have been known for some time for specific models such as the MSSM 
with real parameters UMI] and complex CMS]; and for the NMSSM with real parameters and complex 

[25] , At two loops the momentum effects are expected to be small: according to recent calculations for specific 
models they are comparable to the experimental uncertainty |24H26j and since the momentum-dependent corrections 
due to new physics scale at best as rri 2 H /M^ ew Physics relative to the effective potential contribution, we expect this 
to be a general result. Hence effective potential calculations, with their concomitant great simplification over the 
diagrammatic approach, should be useful at two loop order and beyond (even if the inclusion of the momentum 
dependence will ultimately be necessary to reach the experimental accuracy). 

In general even a two-loop calculation of the dominant contributions at zero external momentum is available for 
just two supersymmetric models: the MSSM |27H43j and partially for the next-to-minimal supersymmetric standard 
model (NMSSM) with real pO] and complex [44] parameters. There continues to be much work in this direction and 
there are now some calculations of the strong (i.e. proportional to a s ) momentum-dependent contributions for the 
MSSM [mf/RlI . These results have variously been made available to the community in model-specific public codes: 
FeynHiggs [55]? SoftSUSY EH. SuSpect [H] and SPheno [9, TD| for the MSSM and NMSPEC m , Next-to-Minimal 
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SOFTSUSY 1181121, and NMSSMCALC [SO] for the NMSSM. There are also some three-loop results, in the Standard Model 
[51], [52] and the MSSM [53] [54] with the code H3m based on [54] . 

The state of the art in these calculations is however somewhat suprising given that the complete generic expressions 
for the two-loop effective potential, valid for a general renormalisable quantum field theory, have been available for 
more than ten years by the work of Martin [35]. These were applied to a complete two-loop calculation of the 
light Higgs mass in the MSSM in the effective potential approach in Ref. [55] ■ Furthermore, generic results for the 
diagrammatic calculation including the momentum dependence up to leading order in gauge couplings have been 
available in the literature for almost as long [561158] , Unfortunately the results of Ref. [39] suffered the so-called 
“Goldstone boson catastrophe” (recently re-explored in PIED]) due to the presence of tachyons in the tree-level 
spectrum so were numerically unstable. Perhaps due to this no public code was made available to exploit these 
prior to Ref. pxP where an implementation in SARAH/SPheno was presented. Currently, the only generic two-loop 
results relevant for the Higgs mass calculation still not present in the literature are the all-electroweak loops and the 
corrections to the Z-boson mass relevant for determining the electroweak expectation value. These will be the subject 
of future work. Here we shall instead continue the process started in Ref. [61] of making the pioneering generic results 
of Martin available in a public code - which entails performing some new calculations. 

As we stated above, calculating the two-loop corrections to the Higgs mass in the effective potential approach is 
expected to be a good approximation. However, there is more than one way to actually perform even this calculation: 
either we can calculate the potential and numerically take the derivatives, as done in Refs. [38] and [61] , or we can 
perform the calculation diagrammatically and set the external momentum to zero. In this work we shall exploit this 
equivalence: we shall analytically take the derivatives of the effective potential, producing expressions equivalent to 
the diagrammatic calculation and having the same structure, but with much simpler loop functions. The advantages 
of this over the first method are that the results are numerically stable 1 ; it is in principle a faster computation for 
more complicated models where the numerical method must make several passes to ensure stability; it can later be 
extended to a full diagrammatic calculation by simply changing the loop functions - but at zero momentum the loop 
functions are much simpler and therefore significantly faster to evaluate. We shall therefore compute the analytic 
expressions for the first and second derivatives of the two-loop effective potential and implement them in SARAH. As 
in Refs. ED El we shall ignore broken gauge groups, and adopt the same ansatze regarding the contribution of the 
electroweak gauge couplings to the tree-level Higgs mass matrix, to which references we refer the reader; the reasons 
for restricting to the so-called “gaugeless limit” are (a) partial circumvention of the Goldstone boson catastrophe 
(complete evasion in the case of the MSSM or any theory where the electroweak gauge couplings entirely determine 
the Higgs quartic potential); (b) significant simplification in the expressions and therefore speed in calculation; (c) the 
electroweak contributions are expected to be small, of the same magnitude as the momentum-dependent contributions. 
On the other hand, in the appendix we provide just the tadpole contributions in the case of broken gauge groups, and 
will return to the full expressions in future work. 

The layout of the paper is as follows. In Sec. [H] we explain our procedure to take the derivatives of the effective 
potential to extract the two-loop tadpole functions for a general theory with massless gauge bosons in a form convenient 
for automation. In subsection |II C| we summarize our results for the tadpole diagrams; we present the full set of 
second derivatives in appendix [B] The implementation of these results in SARAH, including some technical details 


of the translation of the generic results into an an algorithm, is explained in Sec. |III| before we conclude in Sec. IV 


Impatient readers interested in using our implementation of the results might want to jump directly to subsection III B 


II. DERIVATIVES OF THE EFFECTIVE POTENTIAL WITH MASSLESS VECTORS 

In this section we derive the expressions for the two-loop tadpoles in a general quantum field theory with massless 
gauge bosons in Landau gauge. To do this, we analytically take the dervatives of the expressions in [53], Writing the 
couplings in the notation of that paper, the theory is defined by real scalars </>*, Weyl fermions and massless gauge 
bosons A“ where the gauge covariant derivative for the fermions and scalars are 


D^i =8^! + igA^(T a )'jtpj 

Dfttpi =8^ + (Ill) 


i.e. not subject to potential errors from ill-judged step-sizes in the numerical derivation or from parameters being too small. 
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The structure constants are imaginary antisymmetric matrices that obey the gauge algebra but, since we are 
writing in terms of real bosonic fields, for complex representations they will have twice the dimension of the equivalent 
generators T a (so e.g. a U( 1) generator is two-dimensional). We define as usual tr (Q a 8 a ) = )Tb d(i)C(i ) where d(i) is 
the dimension of the representation of field i and C(i) the quadratic casimir, and similarly for T a . 

The lagrangian is then composed of the normal kinetic terms of the scalars and fermions using the above covariant 
derivatives supplemented by purely scalar and scalar-fermion interactions 

C-s = - ljk (f>i(t>j(t>k - 7^ l3kl <l>i<l>j<l>k(t>i 
jC-sf = - l,y IJk 4’i4>j4'k + c.c. (II.2) 

y is in general a dimensionless complex tensor with y IJk = y JIk , while X l ^ k , A* jW are real, symmetric tensors. 


A. Effective potential 



SS 




FFV SSS SSV 


FIG. 1: Two-loop diagrams contributing to the effective potential in the gaugeless limit. 


We can simplify the expression for the effective potential in the Landau gauge given in Ref. 155) for our case; with 
all gauge groups unbroken some diagrams do not contribute. The non-vanishing diagrams are shown in Fig. m 
The contribution of each diagram to the effective potential is given by: 


t/(2) 

V SSS 

= ^(A ljk ) 2 fsss(rrii, m f , m|) 

(11.3) 

t/( 2 ) 

V SS 

=\>? ijj fssimlm]) 

(IL4) 

y( 2 ) 

V FFS 

=\\y IJk ? f FFS (ml, m 2 j , ml) 

(II.5) 

y( 2 ) 

FFS 

= \y IJk y 1 ' J ' k M* n ,Mjj, fppsim], m 2 j , ml) + c.c. 

(IL6) 

y(2) 

V ssv 

= ( jd{i)C(i)fss V {m 2 , m 2 , 0) 

(IL7) 

y(2) 
v FFV 

=Y d ( I ')C( I )f FFV (m 2 I ,m 2 1 , 0) 

(II.8) 

v (2) 

FFV 

= - Y d (I)C(I) m2 lfTD?v( m2 l’ m2 l’°) 

(II.9) 


Here, y and A are the trilinear and quartic couplings of eq. (II.2), g is a gauge coupling, and M 
The loop functions, given in terms of standard basis functions given in appendix [ aJ are the same 


are fermion masses, 
for MS and DR in 
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the case of 


fsss = - I(x,y,z) 
fss =J (x, y) 

Iffs =J (x , y) - J{x, z)-J(y, z) + (x + y- z)I(x, y, z) 
f~FFS = 2-7(x, V} z ) 

fssv =(x + y) 2 + 3(a: + y)I(x, y, 0) +3 J(x,y) - 2 xJ[x) - 2 yJ(y) (11.10) 

but differ for those with vectors and fermions: 



MS 

dr' 

Iffv 

f"FFV 

0 

6 1 (x, y, 0) + 2(x + y) - 4 J(x) - 4 J (y) 

~{x + y ) 2 + 2 xJ{x) + 2 yj(y) 
6I(x,y,0) 


All of these functions are symmetric on the substitution of their first two indices, but may not be so with the third. 
In fact, we can then combine the vector-fermion diagrams to give 

Vffv + V^ v = ^ d(I)C(I)F FV (mj) ( 11 . 12 ) 

where 

Fpy(x) = — 4x 2 + 4 xJ(x) — 6 xl(x, x, 0) + <iMs4xJ(x) (11.13) 

where 6 ms is zero for DR 7 and one for MS. 


B. Derivatives of the potential 


Here we shall analytically take the derivatives of the potential. In a generic model, we may want the derivatives 
of the Higgs potential in terms of some unrotated fields, i.e. in a basis where the mass matrix is not diagonal; let us 
say that we start with such a case. We write the tree-level potential in terms of some expectation values v t and the 
associated real fluctuations Sf as 


V ! 


scalar tree 


= V 0 (t 


0 + 


^ oO qO q0 j_ ^ \ 

g A 0 24 A ° 


ijkl 


qO qO qO qO 
°i D k°l • 


Then there is a tree-level rotation 


(11.14) 


cO _ r>0 q 


(11.15) 


to diagonalise the mass matrix; we then obtain 

V 


scalar tree = y 1 2^ + ^ijkggg + J _ X Hklg g g g 

2 6 24 


(11.16) 


We can write 4>i = v, + Si and work with the couplings of eq. (II.21. We then need the quantities which enter in the 
effective potential calculation, which are masses and couplings depending on {S)}. 

In general, we have 




A ijk {S) 


A ijkl (S) 


=- V 

dSidSj 

=m 2 5 ij + A ^ k S k + ^X ijkl S k S t 

dSidSjdSk 

= X ijk + A ijkl S l 

= yjkl 


(11.17) 
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with the shorthand notation S = {S)}. Hence we can write 

A x “ M(s) =° 


—X ijk (S) =X ljkr (S) 
oS r 

=X ijr + X ijkr S k . 


(11.18) 


Similarly for the fermions we have 


m^j =M 11 ' Mjj, 

=y II ' r M* JI ,+M II 'y JI , 

_d_ 

dS r 


■y Us =o. 


(11.19) 


However, for the purposes of the effective potential, we then require a further diagonalisation for m 2 AS): we put 


Si = R tJ (S) S'. 


( 11 . 20 ) 


We name the couplings in this basis as rhi(S),X z ^ (S),X l:ikl (S). We then express the effective potential by inserting 
the couplings and masses in the basis {S’'} into the formulae of eqs. ( II.3| ) (II.9). However, to take the derivatives we 
rewrite the expressions in terms of the basis {5)}, and use the trick (with m 2 = (m 2 )) 


d 


1 


dS r \q 2 + m 2 


1 

I m 2 


ik 


urn kk' 

dS r 


1 

_l_ m 2 


k'j 


For similar expressions we write by abuse of notation (using C = 167r 2 p, 2e (27r) d [57] ) 


J {m 2 ik ,m 2 jl ) = C 2 J d d qd d k 


ik 


l 


k 2 + m 2 


jl 


For fermion propagators we can write 

M IV 


M. 


1 


1 


m, 


ij 


q 2 + M JK M Kr q 2 + M U M JK 


M, 


KI> 


( 11 . 21 ) 


( 11 . 22 ) 


(11.23) 


Let us demonstrate our method on a brief example: 


1 \ iijj fss (rh 2 , m 2 ) =R° rp ^ r l~X iijj fss(™lrh 2 j ) 

= R rvl^l X ikil fss(m 2 ik ,m§,) 


‘dS r 8 

= R rp\ xlkjl ( S )fss 0) ( m im >™ 2 nk ; m 2 ; ) (m^ n ) 


S->0 


=i^ s A^‘(5)/^(m^,m^;m^)A OTnr (S) 

R r P ]y ksJ X lkr f^f\m 2 , m 2 ; m 2 ). 


(11.24) 


Recall that fss{x,y) = J(x,y) where J is the finite loop function and J(x) = J(x) — f, see Eq. (A.3); here we have 
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defined 


fgf\x, y, z) = — O 2 f d d kd d q-±—-±—^— + - (j (z) Cz f d d k 
J k 2 + x k z + y q z + z e \ J 


k 2 + x k 2 + y J e 2 


= -C 2 

1 


d d k 


1 


1 


k 2 + x k 2 + y 


d d q 


q 2 + z 


—2( z ) + 


x-y 

1 


1 


J(x)-J (y)jj(z) + -J(z) 
J{x) - j(y)]j(z) 


x-y 

= ~ Bq{x, y)J(z). (11.25) 

Note that the R®, a are not functions of Si and so do not present complications if we want to take futher derivatives. 


We can similarly take the derivatives of all the remaining loop functions; these are derived from the basis: 
^-J(m 2 ,m 2 ) -> - B 0 (m 2 ,ml)J(m 2 )^^ - J(m 2 )B 0 (m 2 ,ml)~g^ L 


& 2\ r / 2 2 2 \(d m ik , 


dm 2 . 


^-I(m 2 ,m 2 j,m 2 k ) ->• - U 0 (m 2 , mf, to 2 , - U 0 (m 2 ,mf,m 2 ,m 2 k )-^- 1 


-U 0 ( 


m 2 k ,m^,m 2 ,m 2 )^^ 


and 


d 


dm 2 . 


-7^-m.i I^, m 3 , m k ) , m j ,m k ) - 


dS , 


i drnl 
dS r 

2 

jl 2 tt (^2 _2 _2 _2 


Uo(mf , 


m i i m j , m k ) 


dm 2 , 

dS r 


_2 rj /■ 2 2 „,2 „,2\ kl ' 2 tt 1^,2 ^.2 2 ‘ 2 \ 

dS Wi C/ °^ W ' ’ J ’ Wfe ’ ?7?i ^ , m k , m t , to,-). 


The derivative of a typical term in the effective potential will have the form 

JL A ijk A iikf a{xi , y j:Zk ) =2f a ( Xi , Vj ,Z k )A» k -^ (,!"') 


+ 


/ r)rr? 2 

A V k A i j k ^ f £’ 0 ’°\ Xi , Xi r, yj , Zk) + ( X ^y) + ( x ^ Z ) l 


(11.26) 


(11.27) 


(11.28) 


where, generalising the above, it is straightforward to show that for a generic function appearing in the effective 
potential f a composed of polynomials (even containing monomials with negative exponents) multiplying the loop 
functions above, we can write 


A w >(*, w) s / ‘' (w) - / ‘‘ (w) 


(11.29) 


and similarly for permutation of the indices. On the other hand, this explicit expression is often inconvenient in 
practice due to the need to carefully take the smooth limit when x = u\ it is instead more practical to rewrite 
the right hand side in terms of our basis of loop functions multiplied by suitable polynomials. In the following we 
present explicit expressions for the first derivatives (and, in the appendix, the second derivatives) which have been 
appropriately simplified to remove the apparent singularities. 


C. First derivatives 


Here we gather the set of two-loop tadpole diagrams. There are only three topologies, two of which only apply 
to the all-scalar case. All generic possible diagrams are given in Fig. |II C| Here we present our results for the first 
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FIG. 2: Tadpole diagrams at the two-loop level which don’t vanish in the gaugeless limit. 


derivatives in the basis {Si}, so without the rotation matrices R ; the full tadpole in the original basis is then given by 

m/( 2 ) 

gg 0 = Rr P [Ts + Tssff + T FFF s + Tsv + T fv ] (11.30) 

with the tadpoles on the right-hand side to be defined below. 


1. All scalars 


We start with the purely scalar diagrams which are in the first row of Fig. IIC The entire contribution is given by 


with 


Ts = Tss + T S ss + T S sss 


Tss =-A^ A ikr fgf ] (ml ml ; mj) 

Tsss =Ia"^A ljk fsss(mlml m 2 k ) 

rri 1 \ rii' \ ijk \ i'jk j*(l,0,0) / 2 _2 ^2\ 

Tssss = | A A J A f sss (mi, m t , , to.,-, m k ). 


The new loop functions are defined as 


= - B 0 (x,y)J(z ) 
fsss 0) ( x >Vi u > v ) =U 0 (x,y,u,v). 

Note that f corresponds to Xsss, fsss corresponds to Ssss , and fggg ' 0) to Wssss of Ref. 
zero external momentum. 


(11.31) 

(11.32) 

(11.33) 

(11.34) 


(11.35) 
in the limit of 








2. Scalars and fermions 


We have, first, the diagrams with two scalar propagators: 

Tssff =^y IJk yuifFFs( m F ™ 2 j] m 2 k , ml)X klr 

\y Uk y 1 ' J ' k Mj V M jj,X klr U 0 (ml, ml, mj, m 2 j ) + c.c. 

Then we turn to one scalar and three fermion propagators: 

Tfffs =2Re[y IJr y IKm y KLrn M} L ] T FFFS (ml, m 2 j, m 2 K , rn 2 m ) 

+ 2R e[y IJr y IKrn y JLm M^ L ] T FFFS (m 2 j, m 2 j, m 2 K , m 2 m ) 

- 2R e[y IJr y KLm y MNm Ml K M* JM M* LN \ m 2 , m 2 ,m 2 m ). 

Here we have defined 

fFFs( m i', m *, m J'’ m l ) = - So(mf,,m|)J(m5) + B 0 (ml,,ml)J(m 2 k ) + I(m 2 P ,m 2 j,ml) 

- (ml + m 2 j- m 2 k )U 0 (ml,,ml, m 2 j, m 2 k ), 

/^’](m|,mj;rnfe,mf) =B 0 (ml ,m 2 k )J(ml) + B 0 (ml,ml)J(m 2 ) - /(m 2 ,m 2 ,m 2 ) 

- (ml + to j — ml)U 0 (ml, m\, ml, m 2 ,) 

T FFFs( m /, m J> m L’ m m) =Uo(m F mj, m L , m m ). 

rr ( 2 2 2 2 \ _rl,0,0/ 2 2 2 2 \ 

T FFFS\ m H m Ji m Ki m rn) =JFFS\ m F m Ji m Ki m rn) 

T FFFS (ml,m 2 j,m 2 K ,m 2 m ) =1 (ml, m 2 K , m 2 m ) - mlU 0 (ml,m 2 J ,m 2 K ,m 2 rri ). 


(11.36) 


(11.37) 


(11.38) 


3. Diagrams with vectors 


Finally, the two generic diagrams involving vectors are given by 


T SV =^d(i)C(i)X iir ( 3/(0, m 2 , m 2 ) - J(m 2 ) + 2m 2 


=—d(i)C(i)\ nr ml - 12 + 11 log m 2 /Q 2 - 3log 2 m 2 /Q 2 ) 


=g 2 d(I)C(I)Re(M n ,y II ' r )4m 


6 - 7 log ml/Q 2 + 3 log 2 mj/Q 2 + <%jg[21ogm 2 /(3 2 - l] 


(11.39) 


T fv =g 2 d(I)C(I)Re(M ir y n ' r )4 ^ - 3/(0, m 2 , m 2 ) + 5 J(mj) - 4ml + <^Ms [ 2J ( m /) + m V \) ( IL4 °) 


D. Second derivatives 


To find the second derivatives of the potential we can apply the same technique. However, in principle, we can 
simply use the results of |56l . which computed (diagrammatically) the two-loop scalar self-energies at leading order 
in gauge couplings. Since we want the self-energies for neutral scalars, this comprises all of the contributions, and if 
we want the results in the effective potential approach we can simply set the external momenta to zero (and, for this 
work, the masses of the gauge bosons to zero). In fact, for the majority of the diagrams, these yield the same result. 
However, in a few cases we find that by taking the derivatives of the potential we find simpler results (which are of 
course entirely equivalent). The full result is given by 


dV^ 

dS°dS° 


=RlR" q n i j( 0 ) 


- C>0 7~>0 

- -lX^p q 


n 5 - + n ,SF(W) 


H 


SF 4 (M) 


n S 2 F 3 (M) n S 3 F 2 (V) 

ij ~ r ij 


tjSF^V) 


Hf/+ n r 


(11.41) 
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where the superscripts correspond to the numbers of scalars, fermions and vectors with types of topology listed for 
diagrams (M) and (V). We give the complete set of relevant expressions in appendixjBj the corresponding diagrams 
are shown in Figs. |d| and |dj The expressions for and IIW exhibit particular simplifications in our 

approach. 


/ \ \ 



Ssss Msssss Zssss Ussss 



Yssss Wssss X S ss Vsssss 

FIG. 3: Two-loop self-energy diagrams involving only scalars. 





Mffffs Msfsff MpFFFV 



Vffffs Vsssff Wssff 



Mssssv Wsssv Vssssv 


FIG. 4: Remaining two-loop self-energy diagrams which do not vanish in the gaugeless limit. 
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III. IMPLEMENTATION IN SARAH 

We have implemented the new routines in SARAH. By including the first and second derivatives of the effective 
potential using the analytic expressions here, rather than numerically taking the derivatives of the couplings and 
masses as performed in the previous version |61j . we can guarantee greater numerical stability, accuracy and speed 
improvements. Moreover, this approach allows a straightforward upgrade to the pole mass calculation by simply 
changing the loop functions called to those defined in Ref. [58] based on loop functions implemented in TSIL |63j , 
which will be made possible in a future version. 


A. Method 


For any given supersymmetric 2 model $M0DEL, once the user has specified the particle content and their sym¬ 
metries, SARAH calculates all of the vertices and masses. It then writes a Fortran code (placed in the sugges¬ 
tively named 2LPole_$MQDEL. f90) which implements our expressions, linking to a static Fortran code (named 
2LPoleFunctions. f90) of the basis functions for the generic first and second derivatives of the effective potential 
defined in section |H| and the appendix. These two pieces of Fortran code are called by SPheno during the calculation 
of the loop corrections to the Higgs mass. Here we shall give a few details of how SARAH writes 2LPole_$M0DEL. f 90. 
The overall algorithm is to 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 


Generate masses and couplings for all relevant particles in the gaugeless limit. 

Populate and classify all tadpole topologies according to particle content. 

For each tadpole topology, pass the set of diagrams along with information specifying the symmetries to a 
generic writer function. 


Rotate the total tadpole vector by the Higgs rotation matrix to the non-diagonal basis, c.f. eq. (11.301. 
Populate and classify all second derivative topologies according to particle content. 

For each mass topology, pass the set of diagrams along with information specifying the symmetries to a generic 
writer function. 

Rotate the mass matrix to the non-diagonal basis (as with the tadpoles). 


The writer function is actually identical for tadpoles and mass diagrams with a switch to adjust the number of Higgses. 

It cycles through the list of diagrams and applies the following process for each: 

(a) Determine symmetry factor of diagram due to permutations. 

(b) Determine the colour factor; for diagrams with a gluon propagator this is simply d(I)C(I) whereas otherwise 
we must trace over colour indices of the vertices. In principle, for four-point vertices there can be two colour 
structures for the vertex which superficially leads to more than one colour factor for such diagrams. However, as 
we can simply see by inspecting the expressions in the appendix, or by considering that the colour factors have 
to be inherited from a corresponding vacuum diagram (since differentiating with respect to neutral Higgs fields 
cannot introduce any additional colour factor), for diagrams with a four-point vertex consisting of four coloured 
fields the colour factor is given by a trace over the indices in pairs. Hence such four-point vertices are saved with 
the colour factor of the pairs of indices traced over. To be more explicit, such vertices can only contribute if they 
come from differentiating Vg S which contains the coupling X mn . We then must simply take care that the indices 
of the vertices correspond correctly to the indices that are traced over. 

(c) Write a nested set of loops to sum the diagram over the generations of all particles and, for the inner loop, the 
external Higgs legs, since the most computationally demanding aspect is evaluating the loop functions and this 
can be evaluated before calling the Higgs loop - indeed, we also check that the coupling multiplying it is non-zero 
first too. 


A two-loop Higgs mass calculation in non-SUSY models will become available in future releases 


2 





11 


There are subtleties in translating our results into a form usable by SARAH, both stemming from the fact that in the 
SPheno code the couplings are stored in terms of either real or complex scalars, and four-component spinor fermions, 
while, since our results are based on those of Refs. |55] and [56j and for economy we use real scalars and two-component 
fermions. The translation between the two bases as required by SARAH and SPheno is largely as described in Ref. m, 
however here we have the additional complication for fermions of translating chains of couplings and masses such as 

Ci = R e(y IJr y KLm y MNm M* K M* JM M* LN ). 

In SARAH, the interactions of Weyl spinors ip are derived from the corresponding Dirac spinor interactions U/ as 

6C . ( SC 6C 


Vertex = i- 


Pr = i{c L P L + c R P R ), 


(III.l) 


5(p r 8^i5^j \5(p r 5ipf5ipj 8(p r 8ipj* 5ipj-* 

with ipx = Py^x (V = L, R; X = I,J) and polarization operators Pl , P R ; and so Cl,r O ~y IJr ,c * L1{ O —yij r - 
When we have complex scalars <£>, we should write 

C Dd> m (Tj(c L (/, J,m)P L + c R (I, J,m)P R )*j + $ m (^j (c L {I, J,fn)P L + c R {I, J,m)P R )^j 
= $ m (vj> / ( Ci (j, J, m )P L+CR (I , J : m)P R )^j + $ m (*j(c* R (I,J,m)P L + c* L (I , J,m)P R )*j (III.2) 

and so Cl{I, J, m) = c* R (J, I,m). For a given topology, SARAH populates the diagrams using Dirac propagators (i.e. 
links fermions with its conjugate) and so for the coupling C i above we will find sets of particles 


{*/, *J, <R, Rv, R <T}, {R, R, <P m }- 


(III.3) 


Suppose each of the fermions is a Dirac spinor with Weyl spinors ip R R etc, then to construct coupling C above we 
must sum over the left and right-handed Weyl fermions (which have opposite representations of all gauge groups) and 
thus (noting that SPheno always internally stores the fermion masses as real positive definite) 


Ci = Re 


y^L^HXyi’L _|_ y^R^h^yll>R,1pi,mylt> J R ,ll>^ ,771 




= — Re 

= — Re 


c* R (I, J, r)c* R (N , /, m)c* R (J , N, m) + c L (I, J, r)c L {N , /, m)c L {J, IV, m) 
c R {I , J, r)c R (N, /, to)cr( J, N, m) + c L (I, J, r)c L (N, I, m)c L (J, N , m) 


MjMjMjy 

MjMjMl. (III.4) 


On the other hand, if the fermions are all Majorana then ip\ = ip R and we therefore only have half of this sum, so we 
include an extra factor of 1/2. If we consider another example 


C 2 =Re(y IJr y IKm y KNm M * JN ) 
in SARAH we would generate the set of particles 

{R R, <R, {R, R R}, (R R, Cp m }- 

and thus 


(IIL5) 

(III.6) 


C 2 —J>Re 


y ,mi/ w y y^i R ^Ky m 


Mr 


—> — Re 
= — Re 


c* R {I, J, r)c* L {K , I, m)c* R (J , K, m ) + c L (I , J, r)c R (K, /, m)c L (J 7 K, m) 
c R (I, J, r)c L (K , /, m)c R (J , K, m) + c L (I, J, r)c R (K, /, m)c L (J, K, m) 


Mj 

Mj. 


(HL7) 


The above show that it is straightforward to translate the two-component results into expressions in SARAH 

m n p 

[JI y IiJiSi II l II m iml] 


Re 


L 1=1 1=1 
M 


fc=1 


(—l) m+ "Re 


U Ci(di, Ji, Si) J[[ C R (Ij, Jj , Sj) + 11 C R ( Ij , J / . S' ? j 11 C R (Ij, ./j, Sj) [ Mi k 


1=1 


1=1 


p 


J fc=l 


(III.8) 
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where M = 1 for Majorana fermions and zero otherwise. 

A final point regarding the translation into a basis of real and complex scalars is that the new routines assume 
that there is a unique way of constructing a gauge- and global symmetry-invariant coupling A ?jfc from complex scalars 
other than the complex conjugate of the whole coupling; i.e. if \ ljk is permitted for given complex i,j, k then Ar? is 
not. This is evidently true - if both are permitted then we can generate a holomorphic mass term at one loop which 
violates the premise. However, it is important in that we cannot write gauge singlets as complex scalars if they have 
couplings violating the above condition, no matter how small the couplings - for example for sneutrinos in see-saw 
models. 


B. How to use the new routines 

To study a model with SARAH the general procedure is as follows: the user should download and run SARAH with 
the demanded model. SARAH derives all analytical expressions for mass matrices, vertices, renormalisation group 
equations as well as loop corrections and exports this information into Fortran source code. The Fortran source 
code is compiled together with SPheno and all numerical calculations are then performed by the new SPheno module. 
This includes a calculation of the entire mass spectrum, branching ratios as well as flavour and other precision 
observables [64] . For the mass spectrum all one-loop corrections to any particle are included in a diagrammatic way 
m ■ For a supersymmetric model there are now in addition three options to get two-loop corrections in the Higgs 
sector. The first two are based on the effective potential approach presented in Ref. [61] while the new routines are 
called by the (now default) third option. 

A step-by-step description to obtain a spectrum generator for an arbitrary model $M0DEL implemented in SARAH 
reads as follows: 

1. Download the most recent SARAH and SPheno versions into a directory $PATH. Both packages are located at 
HepForge: 

http://sarah.hepforge.org/ 
http://spheno.hepforge.org/ 

2. Enter the directory and extract both codes 

> cd $PATH 

> tar -xf SARAH-4.5.0.tar.gz 

> tar -xf SPheno-3.3.3.tar.gz 

3. Start Mathematica, load SARAH, run $M0DEL, and generated the SPheno output 

<< $PATH/SARAH-4.5.0/SARAH.m; 

Start ["$M0DEL"] ; 

MakeSPheno [] ; 

4. Leave Mathematica, enter the SPheno directory and create a new sub-directory for your model 

> cd $PATH/SPheno-3.3.0 

> mkdir $M0DEL 

> cp $PATH/SARAH-4.5.0/Output/$M0DEL/EWSB/SPheno/* $M0DEL 

5. Compile SPheno together with the new module 

> make Model=$M0DEL 

After these steps a new binary bin/SPheno$MODEL is available. To run it an input file in the Les Houches format is 
needed. SARAH writes a template for that file which has to be filled with numbers. To enable the new functions for a 
calculation of the two-loop Higgs masses based on our new loop functions the following flags have to be set: 

Block SPhenoInput # 

70 # Skip two loop masses: True/False 

83 # Choose two-loop method 

91 # Gaugeless limit: True/False 
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8 -> 3 chooses the new approach to calculate the loop corrections. The other options for flag 8 would correspond 
the effective potential calculations based on SARAH (8->l for a fully numerical derivation, 8->2 for a semi-analytical 
derivation). Also some hard-coded corrections are available which are based on results in literature: 8->8 uses the 
known ots(ub + a t ) corrections for the MSSM, NMSSM, TMSSM or any variant thereof with up to four neutral CP- 
even Higgs fields and including models with Dirac gauginos [55] , 8->9 uses the corrections of option 8 and adds the 
two-loop MSSM (a t +ab + a T ) 2 results based on Refs. p?5l{5Tl HO] ITT] . Note that the last two options are not included 
by default in the SPheno output of SARAH. To include them, the user must make sure to include in the SPheno.m of 
the considered model 

Use2LoopFromLiterature = True; 

Finally, SPheno is executed by 

> ./bin/SPheno$MODEL $MODEL/LesHouches.in.$MQDEL 
and the output is written to 
SPheno.spc.$M0DEL 


C. Validation 

We have intensively used the SPheno output to validate our new two-loop functions, in particular: 

• We found a numerical agreement of more than 10 digits between our code and using public routines for the 
MSSM based on Refs. [55H571 00; TT for the self-energies. In order to perform this validation, it is necessary 
to use the same assumptions: turn off the first and second generation Yukawa couplings; take the Goldstone 
boson and light Higgs masses in the loops to be zero, and set the tree-level mixing angle of the neutral CP-even 
scalars to a = /? — n/2. 

The excellent agreement between all four possibilites to calculate the two-loop Higgs masses in the CMSSM is 
also shown in Fig. [5jThe parameter point used here is the same as in EQ, 

M 0 = M 1/2 = 1 TeV, A 0 = -2 TeV, tan /3 = 10, sgn(/r) = 1 (III.9) 


• Similarly, we found full agreement with available results for the a s (ab + at) corrections in the NMSSM [20] and 
for Dirac gauginos [531 . 

• We compared the full two-loop corrections, i.e. also including corrections not involving the strong interaction, 
for the NMSSM, models with Dirac gauginos as well as for the B-L-SSM [55] against the results using the other 
two options based on a completely independent implementation in SARAH presented in Ref. m- We usually 
found very good agreement. Tiny differences were based on numerical artefacts in the routines using the effective 
potential ansatz. Similarly, we could reproduce the results of Ref. m for the two-loop contributions to the 
Higgs mass stemming from A-parity violating couplings. 

The new routines of course provide better stability. For example, in the routines based on numerically taking the 
derivatives of the potential, it is necessary to take care with the initial step size; if there are neutral scalars which 
have small expectation values then the results from those methods could become inaccurate - this problem occurs in 
general for any neutral scalar having expectation value Vi <C Msusy- Less significantly, the numerical method can 
suffer (small) errors when there are small couplings present, such that they do not induce a sufficient shift in particle 
masses or couplings upon variation of the Higgs vevs to accurately take the derivative. Hence, it is very important to 
have two independent implementations of generic two-loop Higgs mass calculations in SARAH/SPheno: this is the only 
possibility to cross check results for models beyond the (N)MSSM at the moment. Thus, we highly motivate users to 
test all options for the model under consideration and to compare the results. 


IV. SUMMARY 


We have presented the derivation of a new set of expressions for calculation of the tadpoles and self-energies at 
the two-loop level. These expressions include all generic diagrams which do not vanish in the gaugeless limit and 
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FIG. 5: Comparison between the diagrammatic calculation of the two-loop Higgs masses with SARAH/SPheno presented here 
(diag: diagrammatical) with the effective potential calculation of Ref. [61] (p-rmm: purely numerical, semi: semi-analytical): , 
and with the routines based on Refs. [35H37I |40) [411 (ref: reference). The fixed parameters are those in Eq. (III.91. 


are valid in the limit of zero external momenta. This set of loop functions is simpler than the set of expressions 
obtained by taking the limit p 2 —> 0 in the pole mass functions available in literature so far. This allows for a rapid 
numerical evaluation of the Higgs mass. We have implemented these functions in Fortran and included them in the 
new version of SARAH 4.5.0. This provides the possibility to automatically calculate the Higgs mass in a wide range 
of supersymmetric models with a guaranteed numerical accuracy and stability. The obtained precision for the Higgs 
mass is comparable with the one dedicated spectrum generators provide so far for the MSSM, and can now be applied 
to the study of a wide variety of models. 

Aside from accuracy and stability, one of the principal advantages of this approach is that it is readily extendable. 
It would be straightforward to extend the calculation to non-zero momentum by changing the functions in the code 
and linking with the library TSIL. On the other hand, including the electroweak contributions should be possible by 
applying these techniques to the full effective potential; we presented the expressions in this case for the tadpoles in 
the appendix, but the second derivatives are currently unknown - as are the full set of equivalent expressions in the 
diagrammatic approach. Furthermore, to truly reach the full two-loop precision we would require the two-loop shift 
in the Z-mass that determines the electroweak expectation value. We hope to return to these issues in future. 
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Appendix A: Loop Functions 

The pole masses are constructed from various one- and two-loop functions which are defined in Ref. m- However, 
for the purposes of calculating the effective potential these can be reduced to combinations of standard expressions 
involving just the usual functions I and J. Here we compile this dictionary. We stick closely to the notation of 
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Ref. [57] and make often use of results presented there. Note that we use the standard notation 


In x = log 


x 

Q* 


(A.l) 


where Q is the renormalisation scale. 


a. One loop functions 

At the one-loop level only two functions are needed: 

A(x) =lim[A(ai) + x/e] = x(lna: — 1) 

€—>■0 

B{x, y) =lim[B(ai, y) — 1/e] = — [ dtln[tx + (1 — t)y — t(l — f)s] 
e ^° Jo 

Clearly, we find the following relation to J{x) which is widely used at two loops: (recall C = 167r 2 ^ 2e (27r) _d ) 


A(x) =J{x), A(x) = J(x) = J(x) — 


B(a :,y) —>C f d d k—~ -^ 

J ( k 2 + x){ 


){k 2 + y) y-x 


B{x,y) —>- {J(x) - J{y)) = B 0 (x,y) 

p 2 ^oy — x 


(J(i) - J(j/)) 


(A.2) 


(A.3) 


(A.4) 


where we have introduced the subscript to denote that the external momentum is zero. Let us also denote for future 
use 


Co (x,y,z) = 

which is actually symmetric on all three indices. 


Bq(x,z) - B 0 (x,y) 


y- z 


(A.5) 


b. Two-loop functions 

At the two-loop level we make use of the following set of functions 

J(x,y) =J(x)J(y) 

S 0 {x,y,z) =I(x,y,z ) 
d 

To(x,y,z) = - —I(x,y,z) 

U 0 (x,y,z,u ) = —^—(I(x,z,u) - I{y,z,u)) 
y-x 

d 

V 0 (x, y, z,u) =- —U 0 (x, y, z, u ) 

To{x,y) =lim [T 0 (6,x,y) + B 0 (x,y)ln5\ 

5 —>0 


together with 


M (x,y,z,u,v) — j d d k J d d q 


1 


(fc 2 + x)(q 2 + y){k 2 + z)(q 2 + u)((k — q) 2 + v) 


Mq(x, y, z, u, v) 


I 

(u - y) 

i 

{u - y) 


(U (x,z,y,v) - U (x,z,u,v)) 
(U 0 (x, z, y, v ) - U 0 (x, z, u, v)). 


(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A.10) 

(A-11) 


(A.12) 
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These functions have to following properties: (i) I(x,y,z) is symmetric in all arguments; (ii) To(x,y,z ) is symmetric 
in the last two arguments; (iii) Uo(x,y, z,u) is invariant under the exchange z -n- u and x t— y; (iv) Mo(x t y, z,u,v) 
is invariant under the interchanges (x,z) -H- ( y,u ), (a :,y) -H- (z,u), and (x,y) (■ u,z ). 

An explicit expression for I(x,y,z) is for instance given by [55 


I(x, y, z ) =-(x — y — z)Yny\nz + —(y — x — z)\nxlnz + —(z — x — y)\nxlny 

_ _ _ g i 

+ 2xlnx + 2y\ny + 2zlnz — -(x + y + z) — -£(a;, 2/> z ) 


(A.13) 


with 


£(x, y, z ) =R 2 ln[(z + x — y — R)/2z\ ln[(;z + y — x — R)/2z\ — In — In — 

7T^ 

— 2Li 2 [(2: + x — y — R)/2z\ — 2Li 2 [{z + y - x — -R)/2z] + — 

o 


R =\/x 2 + y 2 + z 2 — 2xy — 2 xz — 2 yz. 


(A.14) 
(A.15) 


Note that x,y < z has been assumed here. 


c. Relations required for the pole functions 


The above loop functions are used as a basis for the various loop functions. However, we also find additional 
combinations such as 


Vsssss(x, y, z, u, v ) =—^—(U 0 {x, y, u, v) - U 0 {x, z, u, v)) 

y-z 

1 1 


( k 2 + x){k 2 + y){k 2 + z) ( q 2 + u)((q + k) 2 + v)' 


(A.16) 


we should compare this to 


V(x, y, z, u) 


( k 2 + x){k 2 + y) 2 ( q 2 + u)((q + k ) 2 + v)' 


(A.17) 


Hence we can write Vsssss(x,y,y,u,v) = —V(x,y,z,u) and remember that Vsssss is symmetric on its first three 
entries. 


d. Simplified loop functions 

For the amplitudes with one or more massless (and possibly identical) fields we find simplified expressions for the 


loop integrals, some of which are collected below. 

1 _ 2 _ _ _ 

I{x,y, 0) =-(-5x -5 y+ {-x + y)ln x + 4ylny + lna;(4a: - 2ylny) - 2(x - y)Li 2 (l - y/x)) (A.18) 

I(x,x, 0) =x(— 5 + 41na; — ln 2 a:) (A.19) 

l _ _ 5 7T 2 

/( x, 0,0)=— -*(lna;) 2 + 2x\nx — -x ——x (A.20) 

Bq(x, x) = — ln.T (A.21) 


It is also sometimes necessary to consider the case of small mass splittings. The results for I(S,x,y), I(8,x,x), 
1(6 , S ', x) in the limits 5 —> 0, S' —> 0 can be found in Ref. [55j and we do not repeat them here. 


Appendix B: Second derivatives of the effective potential 

In this appendix we present the results for the second derivatives of the effective potential. Largely these are 
identical to those in [561 with the external momentum set to zero, but for sake of completeness we repeat the full set 
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here. However, certain expressions become much simpler in this limit, notably some complicated functions involving 
both fermion and scalar propagators, and those involving gauge bosons. 

The full contribution is 


n ij —Ilj j 


n 


SF(W) 

ij 


n 


SF 4 (M) 

ij 


n^(M) + n ^ ( v) 


n 


SFi(V) 

ij 


■n?y + n?y. 


(B.l) 


1. Diagrams with only scalar propagators 


The first contribution, including only scalar propagators, comprises the eight diagrams shown in [3j These are 
unchanged from [56j and are given by 


H fj = \\ ijkl X kmn X lmn Wssss(mlmlm 2 m ,rn 2 n )+ ^ kl X klmm X S ss(m 2 ,rn 2 ,m 2 J 

^_ x ikl x jk my lmnn Yssss{m 2^ ^ ^ ^2 } + X^ mn X klmn Z SSS s(ml W 2 , m 2 m , m 2 ) 

+ lx iklm X^ m S S ss(rn 2 k ,m 2 ,m 2 J + ± (X ikl X kmn + X kl X ikmn ) X lmn \Jssss(m 2 k ,mlm 2 m ,m 2 n ) 

+ ^\^ X inPx™P Ysssss{rn l ml m 2 m ,m 2 n , m 2 p ) + ±\ ikm X* ln X kl *>\ mn *Msssss(m 2 k , m 2 , rn 2 m , m 2 n , rr$B.2) 


Here the loop integral functions are given by: 


W ssss(x,y,z,u) = U 0 (x,y,z,u), (B.3) 

Xsss(x,y,z) = ~J{z)B 0 (x,y), (B.4) 

Y ssss(x,y,z,u) = J(u)C 0 {x,y, z), (B.5) 

Z ssss(x,y,z,u) = B 0 (x,y)B 0 (z,u), (B.6) 

S sss(x,y,z) = - I(x,y,z ), (B.7) 

U ssss(x,y,z,u) = U 0 (x,y,z,u), (B.8) 

Vsssss(x,y,z,u,v) = [U 0 (x,y,u,v)-U 0 {x,z,u,v)\/(y - z), (B.9) 

M sssss{x,V,z,u,v) = -M 0 (x,y,z,u,v). (B.10) 


In the case that y = z, we have the simplification 


Vsssss{x,y,y,u,v) = -V(x,y,u,v). 


(B.ll) 


It should be noted (for example, for the purposes of evaluating the colour factors) that topologies X$ss , Yggg, Z$ss 
arise from differentiating Vg 2 J (and hence the tadpole Tgg) while the others arise from differentiating Vg 2 J s (and hence 
the tadpoles T S ss and T S sss)- 


2. Diagrams with scalar and fermion propagators 

The contributions from diagrams with the topology W are 

n f/ (W ° = \ X z: ’ kl Re[y MNk y M ' N ' l MMM'M NN ']W ss j T (ml,mt,m 2 M ,m 2 N ) 

+-X ljklyMNk y M N i 'W s SFF('rnl, m 2 , m 2 M , m 2 N ), (B.12) 

where we can slightly simplify the loop functions: 

w ssFF(x,y,z, u ) = -2W ssss(x,y,z,u), (B.13) 

W ssFF(x,y,z,u) = -(z + u- y)U 0 (x,y,z,u) - I(x,z,u) + B 0 {x,y)(J(z) + J(u)). (B.14) 
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The contributions from diagrams of the topology M with four fermions are 

n£ F4(M) = Re [y KMl y LNj y K L p y M N p M kki M lu M M m'M nni ] M FFFF S {m 2 K , rn|, m 2 M , m 2 N , m 2 ) 
+2Re[y KMz y L NjyKL'py M Np M LL ' M M m'] M FFFFS {m 2 K , m 2 L ,m 2 M ,m 2 N , m 2 p ) 
+Re[(y KMz y L Nj + y KMj yLNi) yKL' P yMN' P M LL M NN ]M FFFFS {m 2 K ,m 2 L ,m 2 M ,m 2 N ,m 2 p ) 


+2Re \y KMi y LNj yKL P y M N p Mmap Mnn'\ ^- FFFF s{ m %^ to m> to at> m P ) 

+Re[y KMz y LNj yKL P yMN P ] M FFFFS (m 2 K , m 2 L ,m 2 M ,m 2 N , m 2 p ), (B.15) 

where 

M FFFF s( x ^y^z, u ,v) = 2 M 0 (x,y,z,u,v), (B.16) 

M FFFFS (x,y, z,u,v) = (y + z-v)M 0 (x,y,z,u,v)-Uo{x,z,u,v)-Uo(u,y,x,v) + B 0 (x,z)B 0 (y,u), (B.17) 

M FFFFS (x,y,z,u, v ) = {x + z)M 0 (x,y,z,u,v) -U 0 {y,u,z,v) - U 0 (u,y,x,v), (B.18) 

M FFFFS (x,y, z,u,v) = (x + y-v)M 0 (x,y,z,u,v)-U 0 (x,z,u,v)-U 0 (y,u,z,v) + B 0 (x,z)B 0 (y,u), (B.19) 

M ffffs (x , y, z, u, v) = {xu + yz)M 0 (x, y, z, u, v) - xU 0 (z, x, y, v ) - zU 0 {x, z, u, v) 

-uU 0 (y,u,z,v ) - yUo(u, y, x, v) + I(x,u,v) +I(y,z,v). (B.20) 

The results from diagrams of the topology M with three fermions are 

n s 2 F 3 (M) = y km ^ Re [y LN 3y L ' Pk y N ' p ' m M LL ,M NN 'Mpp']M SFSFF (ml,m 2 L , m 2 m , m 2 N , m 2 P ) 

+2Re\y LN -'yLPk‘y N Pm Mjvjv'] M SFSFF {ml, TO i> TO mi to at> m p) 

+Re[y LNj yLPkyNP'mM pp ']M SF sF F (m 2 k ,m 2 L ,m 2 m ,m 2 N ,m 2 P f) + (B.21) 


where 

M SFSFF {x,y,z,u,v) = 2 M 0 (x,y,z,u,v), (B.22) 

M SFSFF {x,y,z,u,v) = (v - x + y)M 0 (x,y, z,u,v) + U 0 (y,u, z,v) - U 0 (x,z, u,v) - B 0 (x, z)B 0 (y,u), (B.23) 
M SFSFF {x,y,z,u,v) = (y + u)M Q (x,y,z,u,v)-Uo(x,z,u,v)-U 0 (z,x,y,v). (B.24) 

The contributions from diagrams of topology V, with three scalars and two fermions are 

nf/ F2(y) = \ ikl \i km (Re[y NPl y N ' p ' m M NN ,Mpp,] V S ssff(™1 rn 2 ,m 2 m , m 2 N , m 2 P ) 

+Re[y NPl y N Pm]VsssFF{ml, m 2 , m 2 m , m%, m 2 P fj , (B.25) 

where 


Vsssff( x ’V’ z ’ u ’ v ) = - 2 Vsssss(x,y,z,u,v), (B.26) 

VsssFF{x,y,z,u,v) = U 0 (x,y,u,v) + (z - u- v)V S ssss{x,y,z,u,v) - (J(u) + J(v))C 0 (x,y,z). (B.27) 


The results from diagrams of topology V with four fermions are 

nf/ 4(V) = 2Re[y KL V < ' M ^y L ' Np y M,N ' p M KK ,M LL/ M MM ,M NN> ]V FFFFS (m 2 K , m|, m 2 M , 
+2R e[(y KLz y K ' M3 + y KLj y K ' Ml )y L N P y M ' Np M KK , M MM ']V FFFFS (m 2 K ,m 2 L , 
+2Rc[y KLl y K ' M:i y LNp yMN' P M KK ,M NN ’]V FFFFS (m 2 K , m 2 L ,m 2 M , m 2 N , m 2 p ) 
+2Rc [y KU y KM jy L Np yM'N P M LL >M MM ] V FFFFS {m 2 K ,m 2 L , m? M , m 2 N , rn^) 
+2R e[(y KLl yxMj + y KLj VKMi) yLN P yM / N , P M MM M NN ] V FF prp S (rn| r ,rn|, 
+2Re\y KLz y K MjyLN P y MNp ]VFFFFs(m 2 K , m 2 L ,m 2 M ,m 2 N , m 2 p ), 


•2 „ 2 \ 
m N ,m p ) 

m M ,m N ,m p ) 


2 2 „ 2 \ 
m M ,m N ,m p ) 

(B.28) 
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where 

V FFFF S (x, y, z, u, v) = - 2V S ssss(x, V, z, u, v), (B.29) 

\ FFFFS {x,y,z,u,v) = - Uq(x, y, u, v) + (v - z - u)V S ssss(x,y, z,u,v) - (J(v) - J(u))C 0 {x,y,z) (B.30) 

Vffffs ( x ! y I ^ u,v) = - 2£/ 0 (x, y, u, v) - 2zVsssss{x, y, z, u, v ) (B.31) 

V FFFFS (x,y,z,u,v) =f FF g 0) (x, y, z; u, v) (B.32) 

V FFFFS {x,y,z,u,v) = - U 0 {x,y,u,v) - U 0 (y,z,u,v) - (x + z)V S ssss(x, y, z, u, v) (B.33) 

V FFFFS (x,y,z,u,v) =f£pf\y,z,u;v) + xf£pg°\x,y,z,u;v). (B.34) 

These represent significant simplifications over the full pole contributions. To recapitulate, 
f FF s°\x,y,U] v) =B 0 (x,y)(J(v) - J(u)) + I{x,u,v) - (y + u - v)U 0 (x,y,u,v) 
fppg(x,y,z,u;v) =C 0 {x,y,z)(J(u) - J(y)) - U 0 (x,z,u,v) - (y + u - v)V S ssss{x, y, z, u, v). (B.35) 


Note that fppg °is symmetric on its first three indices. 

3. Diagrams with one vector propagator 

For self energies of neutral scalars where all gauge groups are unbroken, the diagrams involving one vector propagator 
are particularly simple. 


a. Diagrams with scalars 


We have for diagrams involving scalars 


n 


S J =2 9 2 d(i)C(z) 


X ijkk W S ssv(ml,ml, m 2 , 0) + X ikl X jkl G S s(m 2 k , m?) 


(B.36) 


In [561 the functions are given as (setting the external momentum to zero) 

Wsssv(x, x,x, 0) =31 (x, x, 0) — J(x) + 2x (B.37) 

Gss{x,y) =4yV(x,y,y,0) + 4xV(y,x,x,0) - 2U o (x,y,y,0) - 2U o {y,x,x,0) - 2J(y)B 0 (x,y') - 2 J(x)B 0 (y,x') 
+ 2(x + y)M(x,x F y,y, 0) - 2U o (x,y,y,0) - 2U o (y,x,x,0) + B 0 (x,y) 2 . (B.38) 


However, the expression for Gss greatly simplifies, as we could see by taking the derivative of (11.40): 

Gss(x,y) =2 -Uo{rn%,ml,rnl, 0) - U 0 (m 2 m , m 2 , mf, 0) + B 0 (m 2 , ml) + 2) 


= — 12 


ll(xlnx — y\ny) — 3(xln x — yhx 2 y) 


x-y 


— 2 


Ggs{x, x) = — 1 + 51nx — 31n x . 


(B.39) 

(B.40) 


b. Diagrams with fermions 
For the diagrams involving fermions we obtain 

n£. v = 9 2 d(K)C(K ) [R e{f KL y jKL )G FF (m 2 K , m|) + R e(y iKL yi K ' L ' M KK 'M LL ,)Gpp(m 2 K , m|)]. (B.41) 

Here we have the simpler expressions 

G ff (x, y) =2(x + y)[3U o (x,y,x,0) + 3U o (x,y,y,0) - 5 B 0 {x,y)\ - 61{x,x, 0) - 6I(y,y,0) + 10 J(x) + 10J(y) - 16 {x + y) 
+ <5 Ms 4 [ J ( :r ) + J (y) ~ (x + y)B 0 (x,y)] 


Gppix^y) =4^3[/o(x, y, x, 0) + 3H 0 (x, y, y, 0) - 5B 0 (x,y) - 4 

- ^ ms 4 !! 1 + 2 B o{x,y)\ 


(B.42) 
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Appendix C: Derivatives of effective potential with massive vectors 


In this appendix we shall present the full results for the tadpoles including the possibility of massive gauge bosons. 
To this end, instead of parametrising the gauge interactions via covariant derivatives, we shall instead use the notation 
of fool and supplement our interactions (II.21 with 


£sv = - - \g abi iA»AMi \g abi A^A, 

£fv =9i J A£ip^ I a tl ipj 

with repeated indices summed over, and metric signature (—b++). 


(C.l) 

(C.2) 


1. Effective potential 

The full two-loop effective potential in the Landau gauge was given in Ref. [55] : 

v sss = ^( A l3k ) 2 fsss{rn?,m],ml) (C.3) 

Vs 2 s fssimlm?) (C.4) 

V FFS =\\y IJk \ 2 fFFs(Trij,m 2 j,ml), (C.5) 

Vg s =ly IJ V J ' k M^Mh'fFFs(rnlm 2 j,m 2 k ) + c.c, (C.6) 

V^ v =l(g ai n 2 fssv(mlmlm 2 a ), (C.7) 

V^ s =\g aaU Fvs{mlm 2 ), (C.8) 

V vvs =\(9 abl ) 2 fvvs(ml,ml,m ?), (C.9) 

V FFV =\\9 a i J \ 2 fFFv{rn 2 I ,m 2 j,m 2 a ), (C.10) 

V£> y =\gf gf M 11 'MSj'fpwiml m 2 , m 2 ), (C.ll) 

^gauge =^(5 af ' C ) 2 /gauge(TO 2 ,m^,m 2 ). (C.12) 
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The modified loop functions are 



dr' 

^Mfi 

fssv 

i - A {x,y,z)I(x,y,z) + {x - y) 2 I(0,x,y) 

+ (y - x - z)J{x,z) + {x-y- z)J(y, z) + zJ(x,y) 

+2{x + y - z/3)J(z) 

0 

fvs 

3 J(x,y) 

2 xJ(y) 

fvvs 

4 ^ [(-A (x,y,z) - 12xy)I{x,y,z) 

+{x ~ z) 2 l\ 0 , X, z) + (y- x) 2 I{ 0 , y, z) - z 2 I{ 0 , 0 , z) 

+(z - x - y) J(x, y) + yJ(x, z) + xJ(y, z)] 

+ h J (x) + h.J(y ) 

2 J(z) — x — y — z 

Sffv 

\[(^{x,y,z) - 3z 2 + 3 xz + 3yz))I(x,y,z) - (x - y) 2 I(0,x,y) 
+(x - y - 2 z)J(x,z) + (y - x - 2 z)J(y, z) + 2 zJ(x,y)\ 
+2(—x - y + z/3)J(z) 

-2 xJ{x) - 2 yJ(y) + (x + y) 2 - z 2 

Itfv 

6 1(x,y,z) 

2(x + y + z) - 4J(x) - 4 J(y) 

/gauge 

- 8 x 3 y - 8 x 3 z + 32 x 2 yz + 18 y 2 z 2 )I{x,y, z) 

+(y - z) 2 (y 2 + 10 yz + z 2 )I( 0 , y, z) + x 2 {2yz - x 2 )I( 0 , 0 , x) 
+(x 2 — 9 y 2 — 9 z 2 + 9xy + 9 xz + 14 yz)xJ(y, z) 

+(22 y + 222 — 4:0x/3)xyzJ(x)^ + (x +>• y) + (x z) 

x 2 + 12 yz + 2xJ(x) 

+(x +>• y) + (a: O z) 


where A(:r, y, z) = x 2 + y 2 + z 2 — 2 xy — 2 xz — 2 yz and it is understood that f^ s = f RR + A^g. 


2. Tadpoles 


Clearly our results for the derivatives involving no vectors given in the text are unchanged. However, for all others 
we will have to apply our procedure to the more complicated loop functions and also the derivatives of the scalar 
masses and couplings. In fact, of the new couplings only g abl has a non-trivial derivative, so we require 


and therefore 


m 2 ab (S)S^ = 


d 2 c 


dA-dA- 


m 2 ab (S) = m 2 5 ab + g abi Si + .S, .S' 

g abl (S) = g abl + g abij Sj, 


d 


dS, 
d 

dK niab 


g a bi _ gdbir 


m a b (S) =g abr {S) 


(C.13) 


(C.14) 


In the following we define 


fA’°’°\x,u-,y,z) 


_fa{x,y,z ) - fa{u, y, z) 
x — u 

_fa{x,y,z) - f a (x,y,u) 
z — u 


(C.15) 


i.e. we can give the loop functions used for the tadpoles in terms of those in the effective potential. In general these 
can be simplified, in particular to allow the smooth limit u —> x to be taken, but we postpone that to future work 
where we shall also treat the second derivatives. Here we simply present the full set of tadpole diagrams, modifying 
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(11.301: 


8Vw 
~OSf 


= R°rp[Ts + Tssff + Tfffs + Tssv + Tvs + Tyvs + T FFV + Tpp V + T gauge ] 


where Ts,Tssff,T FFF s are as given in the body of the paper and 


Tssv =^g ai3 g ak3 / 55 y 0) {m 2 , m \; m 2 , m 2 ) + ^g mj g blJ g abr ffsvX m h m p m l: m b) 

Tssv =\ 9 abii g abr f^s 0) (ml,m 2 b -,m^) + \g aaik \ ikr f^\m 2 a -, m 2 , m 2 ) 

Tvvs =\g abl g cb '9 acr fvvP( m li m 

Tffv =2gfg^ J Re(M KF y I ' Ir )f i p^\mj,m%]m 2 j,m 2 a ) + ^gfgl J g abr f FF v\m 2 I ,m 2 J -,m 2 a ,ml) 
t ffv =9i J 9 a /Re{y II ' r Mjj' ) [fey(m?, m 2 j, m 2 a ) + Mjf^^ (rof , m 2 r ; m 2 j , m 2 )] 

+ gfg a /Re(M IK 'M KI 'M*jj,y KKIr )f£f£> (mj, mj ,; m 2 , m 2 a ) 

+ \g a i J gYg abr M n ' Mjj, (mf, m 2 j; m 2 a , m 2 ) 

Tgauge =\g abC 9 dbC g adr f£au£X m l, m d I ^ 2 ). 


1 


(C.16) 

(C.17) 

(C.18) 

(C.19) 

(C.20) 


(C.21) 

(C.22) 
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